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Abstract. The first measurement of the pn — » du> total cross section has been achieved at mean excess 
energies of Q ~ 28 and 57MeV by using a deuterium cluster-jet target. The momentum of the fast deuteron 
was measured in the ANKE spectrometer at COSY- Jiilich and that of the slow "spectator" proton (p ap ) 
from the pd — » p sp duj reaction in a silicon telescope placed close to the target. The cross sections lie above 
those measured for pp —* ppu) but seem to be below theoretical predictions. 
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1 Introduction 



The last few years have seen several measurements of rj 
production in nucleon-nucleon collisions pQ but relatively 
few of ui production 013]. The S- wave amplitude in the rj 
case is strong and the total pp — > ppr\ cross section largely 
follows phase space modified by the pp final state interac- 
tion up to an excess energy Q = yfs — J2f 711 f ~ 60 MeV, 
though there is some evidence for an rjpp final state en- 
hancement at very low Q ^ . Here y/s is the total centre-of- 
mass (cm) energy and m t are the masses of the particles in 
the final state. Quasi-free n production in proton-neutron 
collisions has been measured by detecting the photons 
from n decay and it is found that for Q < 100 MeV the 
cross section ratio R = <Jtot(p n — > pnrj) / crtotipp — ► PPv) ~ 
6.5 Now the dn final state is pure isospin 1 = 0, 
whereas the pprj is a mixture of / = and 1=1. Up 
to Q « 60 MeV the cross section for pn — > dn is larger 
than that for pn — > pnn and this can be understood 
quantitatively in terms of phase space in a largely model- 
independent way jl] • In all meson production reactions it 
is important to have data on the different possible isospin 
combinations in order to constrain theoretical models. It 
is therefore interesting to see whether a similar isospin 
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dependence is found for the w, the next heavier isoscalar 
meson. 

Unlike the n case, the w-meson has a significant width 
(8.4MeV/c 2 ) and so Q is here defined with respect to the 
central mass value of 782.6 MeV/c 2 [7|. The pp — > ppui 
total cross section has been measured at five energies in 
the range 4 < Q < 30 MeV at the SATURNE SPESIII 
spectrometer and at Q = 92 MeV at COSY-TOF g] 
where, in both cases, the u) was identified through the 
missing mass technique. The energy dependence deduced 
is rather similar to that of the rj, except that the phase 
space and pp final state interaction have to be smeared 
over the finite u> width, a feature which becomes important 
close to the nominal threshold 0- 

Attempts to measure the np — * duj reaction using a 
neutron beam are complicated by the intrinsic momen- 
tum spread, which is typically 7% FWHM even for a 
stripped deuteron beam |S]. The alternative is to use a 
deuterium target and effectively measure the momentum 
of the struck neutron. This is made possible by detect- 
ing the very low momentum recoil protons, < 200 MeV/c, 
in the pd — > p sp duj reaction in a silicon telescope placed 
close to the target. Such an approach is feasible at in- 
ternal experiments at storage rings such as CELSIUS or 
COSY because of the thin windowless targets that can be 
used there. Under these conditions the recoil proton can be 
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largely treated as a "spectator" that only enters the reac- 
tion through its modification of the kinematics. The mea- 
surement of the fast deuteron in coincidence would then 
allow us to identify the u> by the missing mass method. By 
varying the angle and momentum of the spectator proton 
it is possible to change the value of Q while keeping the 
beam momentum fixed. The principle of this method has 
been proved at CELSIUS for the pn — > dir° reaction, where 
Q could be determined to 2MeV [§]. 



2 Experimental Set-up 

Our experiment was performed using a deuterium cluster- 
jet target ^0] at the ANKE spectrometer JI] situated in- 
side the COoler SYnchrotron COSY-Jiilich, with the fast 
deuteron being measured in the ANKE Forward Detector 
and the spectator proton in solid state counters. The sil- 
icon telescope used for this purpose is described in detail 
in Ref. |12| and only the principal features will be men- 
tioned here. The three silicon layers indicated in Fig.^ of 
respectively 60 /im, 300 zxm, and 5 mm, covered polar an- 
gles 83° < 8 sp < 104° and ±7° in azimuth. Protons with 
kinetic energies T sp in the range 2 — 6 MeV traversed the 
first layer but were stopped in the second, while those in 
the range 6—30 MeV were stopped rather in the final thick 
layer. Energy resolution of the order of a = 150 keV was 
obtained. The second and third layers were composed of 
strips arranged perpendicular to the beam such that for 
T sp > 8 MeV a resolution of a(8 sp ) < 3° could be achieved. 
For the lower energy protons, neglecting the small non- 
target background, the finite target size led to a(9 sp ) < 5°. 



5mm 200- 




deuterium cluster-jet target 



Fig. 1. Top view of the telescope placed close to the deuterium 
cluster-jet target inside the vacuum chamber of ANKE. Pro- 
tons and deuterons emerging from the beam-target overlap re- 
gion (shaded area) are detected in the subsequent arrangement 
of three silicon detectors (see text). 



There was no difficulty in separating slow deuterons 
from protons via the E — AE method in two ranges: 
2.6 < T sp < 4.4 MeV (70 < p sp < 91 MeV/c) and also 
8 < T sp < 22 MeV (123 < p sp < 204MeV/c). This is illus- 
trated for the lower range in Fig. [21 It is seen here that, 
by choosing the 4.4 MeV upper limit, one avoids the pos- 
sibility of misidentifying deuterons traversing the first two 
layers but missing the third. 
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Fig. 2. Particle identification via their energy deposit in the 
first two detectors of the telescope. Boxes corresponding to ex- 
perimental data are compared to curves predicted for protons 
and deuterons. 



The ability to identify a deuteron in the telescope in co- 
incidence with a proton in the forward detector also allows 
us to obtain simultaneously the luminosity by measuring 
proton-deuteron elastic scattering through a determina- 
tion of the deuteron kinetic energy. For this purpose we 
have calculated the elastic proton-deuteron cross section 
at our energies within the Glauber model |13| . Such an 
estimate agrees with the available experimental pd — ► pd 
data at 2.78GeV/c to within the quoted error of about 
10% The 1% uncertainty in the energy of the recoil 
deuteron, and hence in the momentum transfer, induces 
only a 3% error through the angular variation of the nor- 
malising reaction. Due to uncertainties in the geometrical 
constraints in the target chamber, the acceptance correc- 
tion introduces a 15% systematic error in the absolute 
cross section normalisation. The overall systematic lumi- 
nosity error used to determine absolute cross sections was 
thus taken to be 20%. It should be noted, however, that 
the error in the relative normalisation between different 
beam energies is at most 5%. 

In order to distinguish deuterons with momenta around 
2GeV/c, arising from the pd — > p sp duj reaction, from a 
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proton background that is two orders of magnitude higher, 
inclined Cerenkov counters were installed behind the multi- 
wire proportional chambers and scintillator hodoscope of 
the forward detection system of ANKE 11,15 . To un- 
derstand the detection principle, consider the detector re- 
sponse for a proton and deuteron with the same momen- 
tum. The opening angle of the Cerenkov light cone for the 
faster proton is larger. Thus part of the light can reach the 
photomultiplier after being totally reflected in the counter, 
whereas all the light produced by the deuteron leaves the 
counter. A momentum-dependent threshold was applied 
so as not to change the differential distributions. 

The hodoscope, consisting of two layers of scintillation 
counters, provides an additional criterion for the deuteron 
identification using the energy loss in both layers. By si- 
multaneously varying the AE-cut and Cerenkov efficiency 
level, an optimal combination was found which leads to 
only a 20% loss of deuterons while giving a 92% sup- 
pression of protons due to the Cerenkov counters alone. 
Projecting the energy loss in the second layer along the 
predicted energy loss of deuterons (ex /3J 2 ), one obtains 
the dotted histogram shown in Fig. [31 A further cut on 
the analogous distribution in the first layer reveals a clear 
deuteron peak (solid line) . Moreover, the shape of the re- 
maining proton background can be determined using the 
energy loss distribution of suppressed particles which, af- 
ter scaling, is drawn as the dashed line. This shows that 
the proton background is on the 10% level. 
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Fig. 3. Energy loss in the second layer of the forward scin- 
tillator hodoscope scaled with /3|, where fid is the deuteron 
velocity. The dotted histogram shows the result obtained af- 
ter imposing a 92% proton suppression with the help of the 
Cerenkov counters. The solid lines result from adding a fur- 
ther AE cut using information from the first layer. The resid- 
ual proton background is shown by the dashed line. 



3 Data Analysis 

Having identified a spectator proton in the telescope and a 
deuteron in the forward array and furthermore measured 
their momenta and directions, one can evaluate the miss- 
ing mass mx in the reaction. To clarify the effects of the 
kinematics, it is sufficient to treat the spectator as being 
non-relativistic. To order p^ p we have then 
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E + m d - E d 



(1) 



where fax is the value obtained at p sp = 0. Here p and E 
are the laboratory momentum and total energy of the in- 
cident proton, pd and Ed those of the produced deuteron, 
and rrid and m p the masses of the deuteron and proton re- 
spectively. The square of the pn cm energy, s, can be eval- 
uated purely using measurements in the spectator counter 
and from this Q can be derived: 



s = (m d + m u + Q) 2 w s + 2pp sp cos 9 sp - 



E + m d 



(2) 

where s is the value for a stationary neutron. Because 
the telescope is placed around 8 sp w 90°, ds/d6 sp is then 
maximal and so the value of Q depends sensitively upon 
the determination of the polar angle of the spectator with 
respect to the beam direction. 

Since in our set-up the fast deuteron is measured near 
the forward direction, the same sort of sensitivity is also 
found for mx when using Eq. (1). Now for each beam 
momentum the beam direction could not be established 
to much better than 0.1°, and this may induce a system- 
atic shift of a few MeV/c 2 in the value of mx- On the 
other hand, in view of the to width, the uncertainty in the 
beam momentum (< 1 MeV/c) is unimportant for both Q 
and mx at this level of accuracy. The struck neutron is 
slightly off its mass shell but the off-shellness is controlled 
by the spectator momentum and rests small throughout 
our experiment. 

In Fig. 01 we show our results from the first two sili- 
con layers (70 < p sp < 91MeV/c), where the spectator 
hypothesis should be very good. The angular information 
is important for the missing mass determination but, in 
view of the limited statistics, we had to sum over rather 
wide bins in excess energy. Experience with lo production 
in proton-proton collisions shows that there is consider- 
able multi-pion production under the w peak [5] . Without 
measuring the products of the lo decay, this can only be 
reliably estimated by comparing data above and below 
the lo threshold. Two of the four momenta correspond to 
largely below-threshold measurements and two above, at 
mean values of Q equal to about 28 and 57 MeV. 

There is an indication of a weak lo signal at the highest 
energy and, in order to evaluate its significance, we have to 
master the large multipion background over our range of 
energies. Two different approaches have been undertaken 
to overcome this problem. In the first, pion production is 
modelled within a phase-space Monte Carlo description. 
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The second method is identical to that used in the analy- 
sis of the pp — > ppuj experiment |2] , where the data below 
u> threshold were taken to be representative of the back- 
ground above, being merely shifted kinematically due to 
the changed beam energy such that the upper edges of 
phase space match. This matching of the ends of phase 
space can also be used to check the set-up of the system 
at each momentum. The only significant discrepancy was 
found at 2.807 GeV/c where, in order to account for a 
slight displacement observed in the data, 3MeV/c 2 has 
been subtracted from all mx values at this beam momen- 
tum. As will be shown in the next section, the two different 
analysis methodologies give consistent results within the 
error bars. 

Most of the background can be described by phase 
space convoluted with the ANKE acceptance, which pro- 
vides a severe cut at low mx ■ It should be noted that the 
available np — > dir + ir~ data in our energy range show the 
deuteron distribution to be fairly isotropic in the cm sys- 
tem |16j . In the absence of neutron data, we parameterised 
the total cross section er(s) for the production of N pions 
in proton-proton collisions by 

where sq is the threshold for Nit production. The expo- 
nent pi is fixed by phase space, but A and pi are free 
parameters adjusted to reproduce the pp — > d(Nir) data 
for 2, 3 and 4 pion production |17|. The assumption that 
each of the three contributions follows a (N + l)-body 
phase space, undistorted by A or p resonances, gives a 
description of the mx distributions for different beam en- 
ergies. To model the pn — > d(Nir) background, the energy 
dependence from the pp case has been used to fix the pi, 
with the A being adjusted to reproduce simultaneously 
our experimental distribution at 2.7 GeV/c and the phase- 
space maximum at 2.9GeV/c. The relative normalisation 
between these two momenta was determined from the pd 
elastic scattering data. The adjusted A values, together 
with the relative normalisation established from the lumi- 
nosity measurement, were used to describe the multi-pion 
background at 2.6GeV/c and at 2.8GeV/c, as shown in 
Fig.|3 

Our method gives a plausible description of the back- 
ground under the u> peak at Q ~ 57 MeV but any oj signal 
at Q fv 28 MeV lies close to the maximum of the phase- 
space acceptance and the evaluation of its strength de- 
pends much more critically upon the background assump- 
tions. Nevertheless, within the parametrisation of Eq. (3), 
it is impossible to describe the phase-space maxima simul- 
taneously at the four energies in Fig. 0] without invoking 
some u> signal at Q » 28 MeV. 

To describe the lu contribution to the missing mass 
spectra, we take the pn — > duj matrix element to be con- 
stant over the Q-bin so that the cross section follows phase 
space. This, combined with the decrease of acceptance at 
large Q, means that the mean value of Q is not quite 
at the centre of the bin. Other plausible assumptions, 
such as a constant cross section, would lead to negligible 
changes in the evaluation of the cross section and mean 




m x [GeV/c 2 ] 

Fig. 4. Missing mass mx in the pd — > p sp dX reaction. Two 
of the four beam momenta 2.6, 2.7, 2.807, and 2.9GeV/c 
are largely below and two above the nominal threshold for 
cj-production, with excess energies summed over the ranges 
-58 < Q < -22, -23 < Q < 13, 8 < Q < 44, and 
42 < Q < 78 MeV respectively. The solid line is a multipion 
fit to the data at 2.7GeV/c and 2.9GeV/c, as described in the 
text. There is evidence for an ui signal (simulated as the dashed 
line) at 2.9GeV/c, though the result at 2.807 GeV/c depends 
much more sensitively upon the background simulation. 



value of Q. In the simulation of the pd — > p sp duj reaction, 
the cross section is smeared over the Fermi motion in the 
deuteron using the PLUTO event generator ^8] . This em- 
ploys the Hamada-Johnston wave function |19j though, at 
these small values of spectator momenta, other more re- 
alistic wave functions give indistinguishable results. The 
same event generator is used also for the multipion back- 
ground. 

Turning now to our second approach, the authors of 
ref. noticed that, apart from the cu signal, the shape of 
the pp — > pp X missing mass spectrum varied little with 
beam energy provided that one looked at the distribution 
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with respect to the maximum missing mass. More quanti- 
tatively, if (3 and (i n are cm. velocities at energies T and 
T n respectively, the measured momenta and angles of the 
protons were first transformed, event-by-event, from the 
laboratory to the cm. system with the velocity —f3 and 
then transformed back to the laboratory with the veloc- 
ity +P n - To see to what extent this approach is valid for 
the ANKE spectrometer, which has a much smaller over- 
all acceptance than that of SPESIII we have recon- 
structed the missing mass for the copious proton produc- 
tion pd — » p S ppX. The data at the four different beam mo- 
menta, kinematically shifted to 2.9GeV/c and normalised 
to the same total number of events, are shown in Fig. 
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Fig. 5. Missing- mass spectra of the pd —* p sp pX reaction at 
2.6 (crosses), 2.7 (stars), 2.8 (closed circles), and 2.9 GeV/c 
(open circles) kinematically shifted using the SPESIII proce- 
dure [3]. The data are all normalised to the same total of 100%. 



It is clear from the figure that for mx > 1-4 GeV/c 2 
the shifted data are in mutual agreement at all beam mo- 
menta. For lower missing masses one sees the effect of the 
production of the Z\(1232) isobar, whose position in the 
shifted mass scale depends, of course, upon the beam mo- 
mentum. The figure also nicely illustrates the influence of 
the ANKE acceptance cut, which strongly favours events 
close to the maximum missing mass. 

When the identical analysis procedure is applied to 
the pd — > p S pdX data, the backgrounds away from the u> 
peak at the different momenta are again found to be com- 
pletely consistent. An average background could therefore 
constructed and this is shown for the two above-threshold 
momenta in Fig. [B] The differences between the experi- 
mental data and constructed background shows evidence 
for structure in the uj region and these have been fitted 
to to peaks whose widths were fixed by the Monte Carlo 
simulation. The u> masses obtained from the fits at the 
two momenta, 780 ± 8 and 787 ± 4 MeV/c 2 , do not differ 
significantly from the expected value. 
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Fig. 6. Missing-mass spectra of the pd — > p sp dX reaction. The 
closed circles are the experimental data at 2.8 and 2.9 GeV/c 
whereas the open circles represent the data at the other mo- 
menta shifted using the SPESIII procedure {5J- The differences 
between the data sets (stars) are fitted to the expected oi-peak 
shape to yield the measured production cross section. 



4 Results 

By comparing the residual signal in Fig. 01 with a simu- 
lation of lu production over this range of spectator ener- 
gies and angles, we would conclude from simulated back- 
ground model that <Jtot(pn — > duj) = (2.9 ± 0.8) fib at 
Q = (28±^) MeV and (8.5 ± 2.8) ^b at Q = (57±^) MeV, 
where the uncertainty in Q reflects the total width of the 
bin and only the statistical error in the cross section is 
emoted. The corresponding numbers obtained using the 
SPESIII background technique, (2.2 ± 1.4) yb and (9.4 ± 
3.3) fjb respectively, are consistent with the first method, 
though the statistical errors are larger because we had to 
subtract a background with limited statistics. This con- 
trasts with our first approach where we imposed the con- 
dition that the background should be smooth. Averaging 
the two sets of results, we obtain cr tot = (2.6±1.6±2.3) /ib 
and (9.0 ± 3.2^2 5) A*b at the two excess energies. The sec- 
ond, systematic, error bar includes some contribution aris- 
ing from the ambiguity of the background discussed above 
but others, such as the uncertainty in the luminosity, are 
common to both the signal and background. 

In view of the limited statistics it might be helpful 
to quote upper limits resulting from the fits to the count 
differences shown in Fig.[S] At the 90% confidence level the 
cross sections at 2.8 and 2.9 GeV/c are below (7.5 ±5) /jb 
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and (17 ± 6) respectively, where the second figure is 
the rescaled systematic uncertainty. 

One source of systematic uncertainty comes from the 
restricted angular acceptance of ANKE a problem 
that becomes more serious with increasing Q. The simula- 
tion of the acceptance, illustrated in Fig.EJfor an isotropic 
production distribution, shows that, while the distribution 
is fairly flat at 2.8 GeV/c, few events would be accepted 
close to 90° at 2.9GeV/c. Although at our energies we 
might expect S'-wave production to dominate, when this 
acceptance is weighted with the possible pure P-wave an- 
gular variations of cos 2 9 or sin 2 9, the resulting overall 
acceptance estimate at 2.9 GeV/c is changed by factors of 
of 1.7 and 0.65 respectively These are, however, extreme 
scenarios and a systematic error of half of the difference 
between these values is a generous estimate of this uncer- 
tainty. In more refined experiments, where the statistics 
will allow us to determine the angular distribution, this 
limitation will be avoided. 
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Fig. 7. Predicted angular acceptance for pd — > p sp duj events 
at 2.9GeV/c (solid line) as a function of the deuteron cm. 
angle, assuming an isotropic production process. At 2.8GeV/c 
(dashed line) the distribution becomes more uniform. 
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Fig. 8. Total cross sections for w-production. The pp — » ppu> 
data are taken from SATURNE |2] (open circles) and COSY- 
TOF 31 (open square), whereas our two pn — > dui points are 
given by the closed triangles. Only the systematic errors are 
shown as the statistical errors 1.6 and 3.2 ub at the two energies 
are smaller or comparable in size. The horizontal bars indicate 
the width of the Q ranges. The dot-dashed curve is the semi- 
phenomenological fit given in Ref. |2 to the pp — * ppuj results 
taking the uj width into account. If the ratio for duj to ppuj were 
similar to that for n production one would then obtain the 
solid curve, which predicts a pn — > dui cross section of over 
25 [Jo at 57 MeV. The predictions of the Jiilich group depend 
upon the relative contributions of exchange and production 
current terms and lie between the two dashed curves 122) . The 
only other published estimate f22] is shown by the shaded area. 



compared to 1.53 mb deduced from a compilation of the 
isospin-related pp — > dn + reaction |2"T) . 



The reduction of flux due to the presence of a second 
nucleon in the deuteron target (shadowing) has been esti- 
mated in the ^-production case to be about 5% of the AW 
cross section |2U] and such a correction has been applied 
to our data. These values are shown in Fig. |S1 along with 
those for the pp —> ppuj reaction. 

Due to the momentum distribution in the deuteron, 
the statistics for spectators in the higher range, 8 < T sp < 
22 MeV, are only about a third of those in the lower range. 
Nevertheless, the corresponding missing mass spectra are 
consistent with those shown for the lower spectator ener- 
gies in Fig. with u> cross sections compatible with our 
results in Fig. El 

We have checked our methodology by identifying events 
corresponding to the pd — > p sp dTT° reaction at 1.22 GeV/c. 
Using the same procedures as for the oj analysis, we find 
atot (pn ~* d7r °) = (l-6±0.3)mb at Q = (135 ± 9) MeV, 
where the statistical error is negligible. This value is to be 



5 Conclusions 

In any meson exchange model, the relative strength of 
oj production in pp and pn collisions depends sensitively 
upon the quantum numbers of the exchanged particles. 
If only a single isovector particle, such as the 7r or p, 
were exchanged then, neglecting the differences between 
the initial and final AW interactions, one would expect 
Ctotipn — * pnu>) I a to t (pp — * ppu) — 5. This would explain 
most of the 6.5 factor found in the r\ case 0. Assum- 
ing that the ratio duj to ppu> is as for r\ production, the 
parametrisation of the available pp — > ppuj data [2HS] leads 
to the solid curve, which lies about a factor of three above 
our data. Another estimate is a little lower but similar in 
shape |23) . Both curves lie within the extremes of the pre- 
dictions of the Jiilich theory group [22], where the major 
uncertainty arises from the relative strengths of produc- 
tion and exchange current terms. 
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Taking our 90% C.L. upper limit on the cross section, 17. 
augmented by the corresponding systematic uncertainty, 
would barely bring the data into agreement with the solid 18. 
line of Fig. [SI Even considering only these upper limits, the 19. 
model predictions appear higher than the data. Any the- 20. 
oretical overestimation might be explained if there were 

Or) 

significant isoscalar exchange, perhaps through the u> it- 
self. 

23 

In summary, we have carried out the first measure- 
ment of the pn — > duj reaction by detecting the spectator 
proton from a deuterium target in coincidence with a fast 
deuteron. Although the data are of very limited statistical 
significance, they suggest that the cross section lies below 
the published theoretical predictions. 

In order to clarify the situation further, we are con- 
structing second generation silicon telescopes that will in- 
crease the acceptance significantly. It would then be of 
interest to try to extend this study to the <f> region so that 
one could investigate the OZI rule in the 1 = channel 
to see if the deviations are similar to those in the 1=1 
channel. 
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